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Abstract 

A careful analysis of differential renormalization shows that a distinguished 
choice of renormalization constants allows for a mathematically more fundamental 
interpretation of the scheme. With this set of a priori fixed integration constants 
differential renormalization is most closely related to the theory of generalized func- 
tions. The special properties of this scheme are illustrated by application to the 
toy example of a free massive bosonic theory. Then we apply the scheme to the (/? 4 - 
theory. The two-point function is calculated up to five loops. The renormalization 
group is analyzed, the beta-function and the anomalous dimension are calculated 
up to fourth and fifth order, respectively. 
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1 Introduction 

With the proof of renormalizability of non-Abelian gauge theories in the early 1970's 
the problem of giving a perturbative definition of a renormalizable quantum field theory 
was solved (eg. [|IJ). However explicit calculations in the commonly used dimensional 
regularization are often tedious. This kept the interest in alternative prescriptions alive. 

Quite recently differential renormalization 0, ||, |J has been proposed. For practical 
calculations this renormalization scheme provides two major advantages. Firstly, it allows 
to regularize and renormalize in one step. No explicit regulators or counterterms are 
needed. Secondly, it is possible to keep the spacetime dimension fixed. This is particularly 
useful for dimension-specific theories like the chiral electroweak sector of the standard 
model. 

We start by analyzing differential renormalization from a purely mathematical point 
of view. Differential renormalization is usually formulated in four-dimensional coordinate 
space by writing divergent amplitudes as Laplacians (we restrict ourselves to Euclidean 
signature, x 4 = (x 2 ) 2 , □ = X^<9 2 ) of less divergent expressions. For example, 

1 11, (x 2 \ 1 11, (x 2 \ p 

Ao and Ai are arbitrary integration constants which are kept for dimensional reasons. 

Initially ill-defined integrals are now regularized by the convention that the Laplacian 
should act on the left and the surface term is ignored. According to this rule the singular 
Fourier transforms of x -4 and x -6 can be derived from the well-defined Fourier transform 
of x~ 2 ln(x 2 /A 2 ) (calculated below, Eq. fl55|)), 

f d*xe ix -P 1 f d*x ( ix . A 1 /x 2 \ 1 (p 2 \ 

diff.ren. V U/ V U/ 

r d A xe ix ' v 1 r d A x ( , v „x 1, (x 2 \ 1 2l fp 2 \ 

diff.ren. X 7 X 7 

where A Q n) = 2/(e c A (i)) and (7=0.5772156. . . is the Euler constant. 

The central point in this paper is to exhibit the meaning of the above prescriptions 
for one-dimensional integrals. To this end we perform the convergent angular integrals in 
0) which leaves us with a radial integral / °° dr that diverges at r = 0. We finally 
split this integral into a convergent part J£° which can be evaluated and a singular part 
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Jq 1 which is kept. These entirely well-defined manipulations lead to 

d A x e tx ' p 1 ( f e c \p\ 

4tt 2 x 4 
r/'.r e lx ' p 



4tt 2 x 6 




(4) 
(5) 



Now we compare this result with Eqs. (0) and (|3|) derived by differential renormalization. 
First, notice that the second term on the right hand side of Eq. (^j) has no p 2 -dependence 
at all. To make the right hand side proportional to p 2 we define 



1 dr __ 1 

^3 = ~2 



(6) 



This leads us finally to the equations 



1 - = -i-lnA = -|-lnA 1 . (7) 
o r 2 4 

Eqs. (^|) and (|7|) can be seen as one-dimensional definitions of differential renormalization. 
However Eq. fl7|) naturally relates the renormalization constants via 

In I ^ 3 



For a ratio Ai/A different from exp(-|) the one- dimensional interpretation of differential 
renormalization is not possible. 

We will see in the next section that all ratios of renormalization constants are fixed by 
consistency conditions. Differential renormalization with these a priori fixed ratios will 
be called 'natural renormalization'. 

It was shown |4| that differential renormalization provides a self-consistent definition 
of renormalizable field theories, without referring to the ratio Ai/Ao as given in Eq. (H). 
Moreover in some cases it is convenient to adjust the ratios of renormalization constants 
according to physical requirements ||. In particular for gauge theories it is useful to 
fix some ratios by Ward identities |2], [7]]. However, depending on the gauge, some of 
these ratios may differ from the prescriptions we give. The treatment of gauge theories 
in natural renormalization is still under investigation, first attempts have been successful 

§■ 

The main advantage of allowing for the above one- dimensional reduction and demanding 
Eqs. (^|), (|7]), (§) is that differential renormalization can be understood on a much more 



general footing. We will see in Sec. |2.3| that Eqs. (y) and (|7|) are almost standard in 
the theory of generalized functions. Thus it becomes possible to replace the recipes of 
differential renormalization by mathematically more fundamental definitions. 

In contrast to differential renormalization, natural renormalization is neither connected 
to coordinate nor to momentum space. One has the freedom to choose the most convenient 
representation for the respective problem. 
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The first example where natural renormalization becomes advantageous is the toy 



theory of free massive bosons discussed in Sec. pA[ The mass is treated as two-point 
interaction which leads by power-counting to a non-renormalizable theory in coordinate 
space. With standard differential renormalization it becomes necessary to adjust infinitely 
many constants. It will turn out that these constants coincide with the a priori fixed ratios 
of our approach. This makes it possible to recover the right result immediately within 
natural renormalization. This does not happen accidentally as can be shown in a general 
theorem. 



The main application of this paper will be the y? 4 -theory in Sec. |3.4| . We focus our attention 
to the calculation of the two-point Green's function. It will turn out that the <y9 4 -theory 
performs almost like made for our renormalization scheme: Most Feynman diagrams of 
a given order precisely match into a formula which allows to calculate their sum without 
evaluating single graphs. This enables us to calculate the two-point function up to five 
loops. 

Finally the renormalization group is discussed. The /5-function and the anomalous 
dimension 7 are determined up to fourth and fifth order in the coupling, respectively. 



2 Definition of the renormalization scheme 

2.1 Comparison with differential renormalization 

We start with a generalization of the ideas presented in the introduction. Repeated 
application of the equation 

°/(* 2 ) = i^ 4 ^/(* 2 ) (») 



leads to 



D?l+1 ^ ln (S) = ~ 4n+ln!(n + 1)! ^ forx^O, n = 0,l,.... (10) 

Note that these equations hold strictly only for x 7^ and may be modified by <5(x)-terms 
(cf. Eq. (pO"D). In a renormalizable field theory one needs only a finite number of these 
equations (n = 0, 1 for the <£> 4 -theory) , however it will turn out to be useful to look at the 
general case. 

The function x~ 2n ~ A has no well defined Fourier transform whereas x~ 2 ln(x 2 /A 2 ) has 
(cf. Eq. (|55|)). The differentially renormalized Fourier transform of the right hand side of 
Eq. fllPl) is now determined by the left hand side with the Laplacian translated as — p 2 || , 

d 4 x e ix - p _ P 2n /y\ 
4^^ = pipV^TTj! 11 U 2 J • ( } 

difi.ren. 

We have introduced different renormalization scales A n for each n to stress that they are 
integration constants which a priori are independent from each other and may differ by 
arbitrary positive factors. The A n are interpreted as renormalization scales. 
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Our analysis starts with the introduction of polar coordinates. 

A-K 2 X 2n+A 7T Jo Jo ^ J 

where we have chosen the z-axis to be parallel to p. We evaluate the convergent ^-integral 
and split the r-integral into a convergent part f£° which is evaluated and a part Jq 1 that 
diverges at zero and has to remain unchanged. 

The calculations are in principle straightforward but tedious ||. The result is (Do 1 = 

E? = o) 

d A x e ix ' p " -1 p 2k ( f 1 dr 1 \ 

4tt 2 x 2 ™+ 4 ~ §2(-4fit! (HI)! \Jo r 2n ~ 2k + l + 2n-2k) ^ 

2p 2n (_ r l dr (f\p\ \ _ 1 _ A l\ 

+ (-4) n+1 n! (n + 1)! [~ Jo ~ + R ) ~ 2 (n + 1) ^[k) ' 

From a mathematical point of view we want p-independent integrals to give p-independent 
results. So we are forced to make the following definitions in order to regain the result 



obtained by differential renormalization (IT 



dr 1 



r 



2n-2k+i 2n - 2k 



and (14) 



1 dr 1 n 1 

or 2n + 2 r^i k 

Eq. ( |i~4"| ) is the analytic continuation of the formula 

1 drr n = —^— (16) 
n + 1 

to n < —1. Eq. (15) shows that within our approach we can not equate the renormalization 
scales A n among each other. We find instead 

1 n 1 

lnA n = hiA Y- (17) 

2n + 2 j^Jfe V ; 

for some scale A. The difference InAj — InAj for any i ^ j is a well-defined non-zero 
rational number. If one violates Eq. one changes the definition of convergent integrals 
or generates p-dependences from p-independent divergent integrals. In the differentially 
renormalized y9 4 -theory Ao and Ai are usually equated which however does not destroy 
the self-consistency of the theory since it is incorporated in the freedom of choosing the 
renormalization scheme. 

An overall factor in the renormalization constants is irrelevant, so we choose a renormal- 
ization scale A according to 

"l dr 

— = -lnA. (18) 
r 

Notice that the left hand side of this equation has no explicit A dependence. One 
assumes r~ l to have the implicit local A-term — In A • 5(r) in a similar way as the dif- 
ferentially renormalized version of x~ A acquires the local renormalization dependence 
-\\nA 2 -4ir 2 5^(x) (cf. Eq. (0)). 
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2.2 First results 

'Natural renormalization' corresponds to differential renormalization with the A„ defined 
via Eq. ([171) . It gives a generalization of the usual definition of integrals. 

The renormalization scale A is kept for 'dimensional reasons'. If we integrate over di- 
mensionful parameters then A combines with other ln-terms to provide a scalar argument 
of the logarithms. A is not a cutoff (notice that the integrals over higher order poles fll4|) 




are A- independent), it is neither large nor small (cf. Sec. |3.1| ). 

We summarize the above discussion by giving our definition for the singular Fourier 
transform (A = 2/e c A) 

d A x e tp - x _ p 2n ( fp 2 \ 1 9 v x \ 

i^^TI" (_4) n+1 n!(n+l)! ~ ~ hxV ' 

Moreover we can derive this equation in the spirit of differential renormalization by Fourier 
transforming and translating the Laplacian □ as — p 2 . But then we have to add 5(x)-terms 
in Eqs. (|I|) and (ID) which are now uniquely fixed as 

1 / n 1 1 \ 
-4" +1 "Un+D^+(2E- + — )□"«(*). (20) 

The fundamental divergent integrals (0) and fll8D are easily generalized tcfl 

l (-D m m1 

drr n \n m {r) = / > ' , n ? -1 , m € N , (21) 
o (n + 1) 

/•i , ln m (r) ln m+1 (A) , x 

/ dr ^ = me N . 22 

Jo r m + 1 

So far we have only discussed singularities located at zero. By translation we can shift the 
poles to any point of R. At infinity however one could introduce a new renormalization 
scale Aqo according to 

fOO (jy 

— = lnA co . (23) 
r 

A should be proportional to Aoo for dimensional reasons and it is very convenient! to set 

A = A.. . (24) 

This allows us to generalize Eq. ( ^3|) to 



drr n \n m (r) = - / } — Tl , n ^ -1 , m G N , (25) 
(n + 1) 

- , ln m (r) ln m+1 (A) _ . . 

dr ^ = ^ , me N , 26 

i r m + 1 



1 In fact it is not possible to introduce different renormalization scales A m in Eq. ( p^ ) as can e.g. be 
seen by comparing the (m + l)-fold one-dimensional convolution of |r| _1 with the (m + l)st power of the 
Fourier transform of | r- 1 1 . 

2 By Fourier transforms, e.g., singularities at zero are mapped to singularities at infinity. Eq. ( |l9| ) 
could also be obtained by an (n + 2)-fold convolution of p~ 2 (the Fourier transform of x~ 2 ). In this case 



the integrals are divergent at infinity and our result would depend on A^. Comparison with (19) leads 

to m a. 
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o 



and together with Eqs. (|2ID and (^2|) we get 

drr n ln m (r) =0 . (27) 

All the integrals defined so far can be summarized by the convention 

n = oo n = ^ , InO = lnoo = In A . (28) 

Note that these equations are symmetric under the interchange of zero and infinity which 
comes from the close connection to analytic continuation. 

We close this section with some remarks on changing variables. Integrals that converge 
at infinity may be shifted by definition. However a naive rescaling r i— > ar in Eq. (|I8D 
leads to 

f l / a dar f 1 dr f l / a dr , . A . , , , . . 

/ =/ — + / — = -ln(oA) ^ -In A . 29 

jo ar Jo t Ji r 

To keep Eq. (|18D invariant under rescalings one has to treat the lower limit zero like a 

variable and write Jq 1 dr/r = J V^ dr/r = ln(l/a) — ln(0/a) = — In A. Or, equivalently, 

one rescales the renormalization scale A according to A i— ► A/a. If, like in Eq. ([16|), the 

integral does not depend on A, rescalings do not affect the result. For more complicated 

variable substitutions it is always appropriate to return to the original variables before 

one approaches the limits (cf. the bipyramide graph in Sec. |3.4.4| ). 



2.3 Relation to the theory of generalized functions 

We recognized already in the last section that Eq. flH|) can be understood in the context 
of analytic continuation. In order to include Eq. (Tl8|) into this concept one has to 'care 
for dimensions' and multiply the integrand by the dimensionless factor (r/A) Q , a h- > 0, 

l fr\ a hr 
drr r 



'o \AJ n + a + 1 

which gives (n + 1) _1 for a = 0, n ^ — 1 and a -1 — In A for n = —1. If, according to Eq. 
(p8j)) we replace a -1 by zero we are back at dl8f) . Note that analytic continuation is only 
correct if one uses the factor (r/A) a and if there exists an a- region in C where the integral 
converges. This prescription differs from dimensional regularization by the absence of the 
surface area fl a +i- In general the a-dependence of cannot be compensated by a 

redefinition of the renormalization scale. 

There are other contexts in which we can understand the renormalization scheme. 
Such are contour integrals in the complex plane or lattice theory which generalizes the 
Riemann sum prescription and eventually provides a purely numerical definition of diver- 
gent integrals 0. Here we present the relation to the theory of generalized functions. 

Assume we are interested in an integral which contains the generalized function / that 
is given as derivative of another generalized function F' = f . With a test function if we 
obtain (e.g. |§) 

/+oo /"+00 
dxf (x) V (x) = (/, If) EE (F\ If) EE (F, if') = - dxF (x) If' (x) . (30) 

-oo J — oo 
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If (p is sufficiently constant at the poles of F the right hand side converges and can be 
used to define the integral on the left hand side. 

Let us e.g. take f(x) = x A 0(x); Q(x) = 1 for x > and Q(x) = for x < 0. We 
choose (p = 0(1 — \x\) where the edges at x — ±1 may be smoothed to be C°°. In the 
limit where this becomes irrelevant we have for A / -1 

pi poo ™A+1 j 

(/» = / dxx x = - / dx— — (-S (x - 1)) = — — (31) 



/o Jo A + l v y " A + 1 

which coincides with our renormalization rule (|16|) . For A = — 1 we may take F(x) = 
ln(x)6(x) and get 

1 Arp roo 

dxln(x)(-8(x-l)) =0 (32) 



OX JO 



which is Eq. ([18]) for A = 1. The same holds for Eqs. ( pi]) and (|22|). 

To see what happened with the renormalization scale A we have to notice that the 
above calculation is ambiguous. There exist several functions F which have the same 
derivative. On the real line they differ by a constant which is irrelevant since the test 
function ip vanishes at ±oo. 

In general however the number of undetermined parameters equals the number of 
disconnected pieces of the integration domain. A singularity of the integrand / at the 
origin splits R into two disconnected parts R~ and R + . Each of the functions F(x) + 
C + DQ(x) is with the same right an integral of f(x) on the real line. However they give 
different results to the integrals ([H]) and fl32|), 

Jdxx x = -J dx lc x + (x) f ^— + D x j J {5(x + l)-5{x- 1)) = j^-j + D x , (33) 



A + 1 JJ A + 1 

J— = - jdx (C_i + 6 (x) (lnx + D_ x )) (5 (s + 1) - <J (x - 1)) = Z^i , (34) 



o 

i 



where all the Ca and D\ can be chosen separately. 

The way out of this ambiguity is to change the topology of the integration domain. 
We can compactify R to R by adding oo (= — oo). Since for A < — 1 the integral over x x 
is well-defined at infinity we can define F over R\{0} which is again a connected domain 
with one integration constant. The D\Q(x) term is discontinuous at infinity and thus no 
longer allowed in Eq. (|33|) . The integral Jq 1 dxx x acquires again the unique value (A + 1) -1 . 

For A > — 1 the integral diverges at infinity and the gluing is not possible. However 
for A > — 1 the integral is finite at x = and can be defined on the connected domain 
R. Just for the case A = — 1 the ambiguity remains since the integral is divergent both at 
zero and infinity. The function In \x\ can only be defined on R~ UR + and one should keep 
the arbitrary constant 7J_i in Eq. fl34|) . With the more intuitive relabeling 7J_i = — In A 
we are back at Eq. fllsp . 

In practice the introduction of the renormalization scale A is a matter of convenience. 
It will turn out to be useful to have this parameter at hand. In principle one could set 
A = 1 using the standard theory of generalized functions and recover A in physical results 
by getting the dimensions right in logarithmic terms. 
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3 Applications 

Now we turn to physical applications. In the following we are mainly concerned with 
four- dimensional integrals which we normalize according to 

/ dx = [ J^L and 5(x) = (2tt) 2 tfW (x) . (35) 

This eliminates all irrelevant factors of 7r from the theory. We get e.g. / dxS(x) = 1 and 
/ dxe tpx = 5(p). Analogously n-dimensional integrals are normalized by (2ir)(~ n / 2 \ The 
metric is always Euclidean. 



3.1 A toy example: the free massive bosonic theory 

As a first test let us calculate the four- dimensional free massive boson propagator in 
coordinate space. The result is well known, 

/Q — Vp-X 
dp^> 9 = TT-^l i m \A) 
pA _|_ m A | x | 

The propagator A is perfectly well-defined. However it is not analytic at m = since the 
series contains logarithmic terms in m. 

Now let us treat the mass (— m 2 ) as a two-point interaction and study perturbation 
theory around m = 0. 

• • = • • + • *- 2 • + • *- % *- 2 • + • • • (37) 



The free propagator in four dimensions is 



dp\e~ ipx = \, (3t 



p 2 x 2 



and therefore 



A (x) = — — m 2 [ dxi 5-— o + (—m 2 ) [ [ dx\dx 2 n 9-^7 + . . . . 

K ' x 2 J \ X - Xl f x 2 V JJJ 1 2 { X - Xl f ( Xl -x 2 ) 2 x 2 

(39) 

Since the 'coupling' has mass-dimension the terms become more divergent with every order 
and the expansion is non-renormalizable. However we can treat the integrals according 
to our rules and obtain an unambiguous result which contains by construction only one 
renormalization scale A. 

In order to evaluate the integrals we can use D(x — x\)~ 2 = — 8(x — xi) and Eq. (|]) 
to derive a recursive formula. Here it is even simpler to remember that the n-th term is 
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the Fourier transform of p 2n . Eq. ( Jl9| ) gives 

A na t. rC n. {x) = ^ + (-W^)^ f dp 



p2k+4 



e %p-x 



k=0 

1 b'"(b¥«)'/ (W\ t „ , t „ A ,„„, 

= ^ + -£Wrirl~)"* ( ' c+2, "* (fc+1) )- (40) 

One could study the renormalization group by looking at rescalings of A. In fact com- 
parison with Eq. (|36|) shows that the situation is even simpler. We just have to equate 
A = m to obtain precisely the correct result. This does not happen accidentally as we 
will see in the next section. 

The differentially renormalized result can be obtained by using Eq. (|TTJ) instead of Eq. 
(0), 

ff ' ren ' = ? + t £ Wriy r H~ J " 2C j • (41) 

It is necessary to adjust the infinitely many parameters A& precisely according to the a 
priori settings ( |TT| ) of our scheme. 

Dimensional regularization leads to a series in 4 — n with a simple pole, 

1 m ! " (mV/4)'/ 2 , /AV\ t „ A „„ 

Since the series is non-renormalizable, it is not possible to renormalize by introducing a 
finite number of counter terms. If one nevertheless tries to follow a minimum subtraction 
prescription, one misses a term — + to obtain the correct result. In the next section 
we will present a general method that allows us to calculate this term. 



3.2 A theorem on singular expansions 

Let us summarize what we did in the last section. We started from a well-defined integral 
/ dpexp(—ip ■ x)/(p 2 + m 2 ) which we tried to expand into a series at m = 0. To this end 
we expanded the integrand into a power series / dpexp(—ip ■ x) J2kLo(~ m 2 ) k p~ 2k ~ 2 . The 
interchange of the sum and the integral led us to the perturbation series in coordinate 
space J2T=o(~ m2 ) k I dpexp(—ip ■ x)p~ 2k ~ 2 . This interchange is obviously illegal. Firstly, 
the integrals diverge at p = 0. Secondly, we obtain a power series in m and we know that 
the correct result has no such representation but contains logarithmic terms in m (Eq. 
(|36l)). Although the integrand is analytic at m = the integral is not. So necessarily the 
expansion is wrong and the diverging integrals reflect this fact. We want to study the 
issue how to reconstruct the true result from such an incorrect, singular expansion. 

Let us slightly generalize the situation and look for the expansion of an integral 1(a) = 
J dxf(x, a) into a series at a = 0. The integrand has a Taylor (or Laurant) series f(x, a) = 
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J2k a k f( k '{x, 0)//c!, but in general we can not expect that the series of 1(a) is given by the 
integrals over the coefficients f( k '(x,0) since the integrals may diverge. We define 

k 

AI (a) = J dxf (x,a)-E^|/ dx f {k) ( x , °) > ( 4 3) 

and conclude that AI(a) will only be zero if 1(a) is analytic at a = 0. So A I (a) gives the 
part of the expansion of 1(a) that can not be reached by standard perturbation theory. 

We call AI the non-perturbative part of the expansion. A priori we know almost 
nothing about it. However in many cases where 1(a) is not analytic at a = one can 
calculate AI (a) by the following theorem. 

Theorem. Assume the integrals J dxf( k \x,0) are regular at x ^ 0. If there exists 
a neighborhood of x = 0, a = where / can be written as / = J2e fi with fi(x, a) 
integrable at x = and the f# k \x, 0) having the following properties 

ft ] (\x\, 0) oc |x| n(M) ln m(M) (|x|) with (44) 
given I : lim n(k,£) jk < or (x, 0) = for almost alii k and (45) 

fe— »oo 

given k : n (k, I) > for almost all I , (46) 

then 

AJ(a)=X; Jdxf e (x,a) . (47) 



Note that all the integrals may diverge and have to be defined according to the rules given 
above. The range of integration can be M n or M + , subsets can be taken into account by 
using step-functions. 

Proof (sketch). Without restriction we can assume that the support of / is a little ball 
B £ around x = since the integrals over the remaining domain are regular and therefore 
do not contribute to AI. Moreover we can assume \a\ to be small, so that we can write 
/ as a sum over fe. Since the fe are integrable at x = one gets for sufficiently small e 

/ dxf (x, a) = dxf e (x,a) 

l JR" ^ JR"\B E ^ k\ 

In the second integral the singularity at x = is excluded and fl4"5|) assures that the sum 
can be interchanged with the integral for small enough a, yielding 

i **i (*. -) = E I */« (*. -) + E E a i (i - U ^ («. o) • 

Now we can use the central argument of the proof. The last integral over the entire R™ 
vanishes since (Eq. fl44|)) it is proportional to dr r N ln M (r) and all those integrals are 



3 all up to a finite number 
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zero in our renormalization scheme (Eq. (0))- We finally use Eq. (|46| ) to interchange the 
second sum over £ with the sum over k and the integral. Therefore 



AI (a) = 



[ dxf (x,a)-J2jr I dxf ik) (x, 0) = £ I dx h ( x > a ) 



□ 

Now let us use the theorem to derive AI(m 2 ) of the scalar bosonic theory. We get 
f(p,m 2 ) = e~ w ' x /(p 2 + m 2 ). Expanding the exponential yields f(p,m 2 ) = J2e{~W ' 
xY/(£\(p 2 + m 2 )) = J2efi(p,m 2 )- Moreover /j fc) (p,0) = (-ip ■ x)\-\) k k\p~ 2 ~ 2k j £\ and 
f$ k \\p\,0) oc \p\ i - 2 - 2k meets Eqs. (M), (M), We can apply the theorem and obtain 



Ar / 2 n ~ rd 4 x {-ip-xf ^(-i|a;|)'l f . . t [°° dpp e+3 

(4 



, n ., Air 2 £\ (p 2 + m 2 ) £\ ir Jo Jo p 2 + m 



2 ' 



The ^-integral vanishes for odd £. The divergent p-integral can be reduced to fundamen- 
tal integrals as follows: j °° dpp 2e+3 / (p 2 + m 2 ) = j Q °° pdp((p 2 + m 2 ) — m 2 ) l+l /(p 2 + m 2 ). 
With Eq. (|27|) we obtain (-m 2 Y +1 J °° pdpj \p 2 + m 2 ) = ±(-m 2 ) £+1 ln(p 2 + m 2 )|g° = 
|(— m 2 )^ +1 ln(A 2 /m 2 ) (with Eq. (p8|) a). The result is proportional to ln(A/m) and van- 
ishes thus for A = m. This confirms the explicit calculation of the last section. 

The above theorem holds for any spacetime dimension. Hence it should as well be 



possible to apply it to the dimensionally regularized result and 'correct' Eq. (42) by 
adding the non-perturbative part. We start with the n-dimensional analogon of Eq. 
All integrals are standard and one obtains 



AT ( A HNfjtl . 4 -n fin ■ n-2 q I q f* ' <W 



-n—l 



£\ (2vr)" /2 Jo p 2 + m 

m 2 (- x 2 m 2 ) e frn 2 ^ {n ' m /4 



£=0 
oo 

r ^ I I L I Oj IIV I I lib I 

2 l + 2 £\ [zK 2 ) ' 

~ m 2 ^+V £ / 2_ 

f^^+H\ (£+1)1 V 4-n 




Together with A dim reg . (Eq. (|4"2"D) the renormalization scale drops out and one obtains 
the full propagator. 

So in the example of a free massive bosonic theory we do not have to go through the 
standard renormalization business. One can use the above theorem instead. The simplest 
way to expand the propagator is using natural renormalization, however dimensional 
regularization leads eventually to the same result. 

In a realistic field theory with dimensionless coupling the situation is slightly different. 
The path integral is a priori ill-defined and the renormalization scale an intrinsic parameter 
of the theory (like e.g. in the integral Jq — = — In A). It makes no sense to equate 



4 More precisely / °° pdp/(p 2 + m 2 ) = f Q pdp/(p 2 + m 2 ) + dp(p/(p 2 + m 2 ) — 1/p) + dp/p = 
±ln(l + to -2 ) - |ln(l + to 2 ) +I11A = InA/m. 
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the renormalization parameter A with the coupling. However it is challenging to try 
to generalize the theorem to path integrals providing a non-perturbative but analytic 
definition of a quantum field theory. 

Anyway the theorem on its own is useful in many elementary mathematical applications. 
Integrals like jj dx/(x n + a"), jj dx\n m (x/A)/(x + a), ^ dxe~ bx / (x + a) n , / °° dxe~ bx - a / x } 
etc. can be expanded at a = by virtue of the theorem ||. 



3.3 Fourier transforms 

Before we start to study <£> 4 -theory it is useful to discuss Fourier transforms since many 
Feynman amplitudes are determined by multiplications and convolutions. 

To this end we generalize the Fourier transforms discussed in the beginning (Eq. flT^)). 
It is convenient to derive the result by analytic continuation. A straightforward calculation 
gives (A = 2/(e c A)) 



d 4 x e ip ' x ( x 2 \ a p 2n ( p 2 \ a e 2Ca T (— n — a) __ p 2n ( p 2 \ a 



(2tt) 2 x 2n+i \A 2 ) 2 2n + 2 \A 2 J T(n + 2 + a) 2 2n+2 \A 2 / 



a n/ a e , (50) 



~ r(-n-a) r(l + a) /~ 2( {2k + 1) 2k+1 \ 

where > a nl a l = — ^ ^^rr^ ^exp > -a 2h+L . (51) 

^ n * r(l-a) T{n + 2 + a) 1 \^ 2k + 1 J K J 

Now it is easy to determine all Fourier transforms of the form x~ 2n ~ 4 ln m (x 2 /A 2 ). To 
produce the logarithms we divide Eq. ( |50|) by a m and pick up the finite term in the 
a-expansion: 

ln m (x 2 /A 2 ) ip „ p 2 " - /y\ 

dx x 2n+4 e =(-!) m! ^LVft-m^j> (52) 

where gi(p 2 /A 2 ) is the finite term of (p 2 /A 2 ) a a l at a = 0. We obtain ^ = ln'^(p 2 /A 2 )/|£|! 
for I < 0. For £ positive gi(p 2 /A 2 ) = for all p 7^ 0. A more detailed calculation || shows 
that p- A gi (p 2 /A 2 ) = \8{p). 

Let us take e.g. n = 1 yielding E^-i a u a l = \a~ Y - 1+ fa+ (-§ + |C(3))« 2 + C(a 3 ) 
and therefore (cf. Eq. (|19l)) 

dse**^ = -p 2 f - In f t-\ - -) , (53) 
x 6 16 y \2 \A 2 J 4/ ' V ; 

, , ln(x 2 /A 2 ) 1 2 / l l2 ^ 2 \ 5, /p 2 \ 17\ 

^•■-4^ = Tq p ("4 (a 2 ) 4 (a 2 ) " ~8J • (54) 



n 



-1 gives ££Li a_ M c/ = 1 + |C(3)a 3 + 0{a 5 ), thus 



/ 



^,4 1 , / fe ,->(*y A2 > = X^f) . ( 55) 



x 2 p 2 J x 2 p 2 



With n = —2 we finally obtain the standard formula / dxe ip x = 5(p). Less obvious is 
/ dxe ijhX ln(x 2 /A 2 ) = -4p~ 4 + 5{p). 
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(a) 



a 



(b) 



(c) 




(e) 



-@ 

(/) («) 



(A) 



(0 



(it) 



(/) 





i r 3 



1 2 




(«) 



^ 

0>) (<?) 



^ — 0^ ^ 

(r) (s) (t) (u) 

Figure 1: Feynman diagrams of the yj 4 -theory. 

3.4 The massless ip 4 -theory 

The first serious test of the renormalization scheme is the discussion of the y) 4 -theory. 
Note that once the Feynman rules and the propagator are fixed the results are unique. 
There is no freedom to choose a certain subtraction scheme. 

We keep our integral normalization of {2n)~ 2 which results in a rescaling of the coupling 
by (27r) 2 . So g is related to the usual 'irrationalized' coupling via 



9 
4 



A 



16tt 2 



(56) 



The Feynman diagrams we are concerned with are depicted in Fig. 1 (a), . . . , (u). The 
corresponding amplitudes are labeled by G a , . . . , G u . With natural renormalization we 
have the freedom to switch between coordinate and momentum space. However, most 
often it is convenient to start the calculation in coordinate space where, at least at higher 
loops, the Feynman rules are more transparent. The final result is given in momentum 
space to make it easier to compare it with other work. 
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3.4.1 Simple results 

The free propagator (a) is given by G a = p~ 2 . 

The loop in the diagram (b) gives rise to a term / dpp~ 2 in momentum space or a term 
(x — x)~ 2 = 0~ 2 in coordinate space. Both expressions are set to zero in our scheme: 
Gb = 0. In this aspect it behaves like dimensional regularization. 

More generally, diagrams that contain tadpole insertions give zero and can be dropped. 
This remains true for any number of internal lines the tadpole may have: G c = 0. The 
reason is that in a massless theory a tadpole insertion can only give rise to a number times 
a momentum conserving (^-function. On the other hand it has dimension p 2 and the only 
number with this scaling property is zero. In our renormalization scheme A 2 occurs only 
in combination with logarithms. 

Moreover, due to translation invariance and Eq. (Wn ), all vacuum bubbles vanish: 
G u = 0. So only connected diagrams contribute to the two-point function. Altogether 
this reduces the number of relevant Feynman diagrams considerably. 

Diagram (d) for n = 1 is the sunset diagram. It was already calculated in the last sec- 
tion. The triple line gives x~ 6 which transforms into momentum space as (|) 2 (| ln(p/A) 2 — 
|). Together with the two external legs and the symmetry factor | we obtain 



3.4.2 Chain graphs 

We call diagrams of type (d) chain graphs. To any order there exists one chain graph 
and, if we disregard the vanishing diagrams with tadpoles, the only remaining diagrams 
up to three loops are chain graphs. 

It is possible to calculate chain graph amplitudes for any n by Fourier transformation. 
In momentum space the series of bubbles gives rise to the n-th power of ~(— ln(»/A) 2 + 1) 
(cf. Eq. (|Hf)). A final convolution with p~ 2 (use Eq. @ as suggested in Sec. |3TT| ) provides 
the result as an n-th order polynomial in ln(p/A) 2 with purely rational coefficients. In- 
cluding the symmetry factors and the external legs we obtain for n > 2 (the case n — 1 
has an extra symmetry which changes the symmetry factor from 2~ n to |) 



There is a nice way to compile this result by a generating function. If we multiply G^ n 
with the factor a n = 3 n ~ 2 for n > 2, a\ — 1, it reproduces the leading logarithms of the 
full </9 4 two-point function correctly. The result may be seen as some approximation to 
the propagator. We get 




(57) 




(58) 




(59) 
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We easily read off 

45 
32 





G dA = [±] ^|^ln 4 (^]-^ln 3 (^U^ln 2 |^|-^li.|^ 




3.4.3 Four loops 

Before we start with the analysis of four and five loops a word of caution is in order. In 
general it is not sufficient to define the integral over generalized functions for defining a 
field theory, since also products of generalized functions appear. In principle e.g. one has 
to consider terms like x 2 5(x) since they might give finite contributions after multiplication 
with x~ 2 . 

In the following we do not care about such terms. The main message of the next 
two subsections is to show that there is a miraculous matching of Feynman amplitudes 
in the natural renormalization scheme that makes calculations easy. This matching is 
not affected by the above problems nor are the leading logarithms of the results. This is 
confirmed by the existence of the renormalization group equation studied in Sec. |3.4.6| . 



Diagrams (e), (/), (g), and (h) remain to be evaluated. G e is basically the square of G^x- 

Gf can be calculated with Fourier transforms. Adding propagators from the interior loop 
to the exterior lines we obtain (see Eqs. d53|) -(|55D) 



i ^ pVii (py 5 \ p- 4 1 f 1 ] (py_s\ t i ( i, i^_y_^ 

x * ' 4 2 1 2 n VA7 4 J UVUAA/ 4 J ' 4 2 x 2 I 2 n \AJ 4 







> 












4 2 x 6 \ 2 \AJ 4/ 4 4 \8 \AJ 8 \AJ 16 8 \AJ 16 
Together with the external lines and the symmetry factor one gets 



/ W p 2 \96 \A 2 J 48 V A 2 ; 32/ v ; 

We are left with two diagrams each of which cannot be calculated by Fourier transfor- 
mation. Both have the same symmetry factor \. This makes it possible to use a formula 
which is specific to four dimensions. 

3 + m - a ♦ \w (65) 

X\ X X\ X Xi X Xi 
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This equation holds up to a total derivative proportional to 

d ( 1 x>*-x% 1 x^-x^ 1 



dx 11 \ (x — X\) (x — x 2 ) (x — x 3 ) (x — Xi) (x — x 3 ) 



2 



(66) 



(notice that gfjr j^z^kz = \8{x—Xx)). After integration over x the total derivative vanishes 
and we obtain: 



+ ■ < ft > . = + ±-^=„ l =^ (67) 

Xi\\Js %2 Xy X X$ %2 2 X\ X3 %i 

Thus we need not solve each of the complicated diagrams (g) and (h) separately!. Their 
sum is equal to two chain diagrams. 

Eq. (|55|) can be interpreted as integration by parts which also proved to be useful 
within dimensional regularization [0] . However only in natural renormalization it allows 
one to calculate the sum of diagrams without evaluating single graphs. One should take 
this as a hint that calculating single diagrams is in general not an appropriate method to 
evaluate higher order perturbation theory. All diagrams (or at least groups of diagrams) 
of a given order should be treated as a unit and calculated together. This strategy will 
be even more useful in the next section. 

3.4.4 Five loops. 

Apart from the trivial diagram 



4/ p 2 V 96 V A 2 J 64 V A 2 I 24 \A 2 J 128 



and the five loop chain graph, Eq. (0), ten diagrams have to be evaluated. These graphs 
split into three classes: (1) Diagrams that can be solved by Fourier transforms (j - I). 
Let us call such diagrams Fourier graphs. (2) Diagrams that can be reduced to Fourier 
graphs via integration by parts (m - r), and (3) the nonplanar diagram (s) that we call 
the bipyramide graph. 



Fourier graphs. Every graph that reduces under the replacement of multiple lines 
(•=•, • w ) and iterated lines (• — • — •) by a simple line (• — •) to the free propagator 
can be solved with Fourier transforms. The calculations are analogous to the evaluation 

5 We do this in the next subsection. The single results will be more complicated than the sum. Each 
of the amplitudes (g) and (h) has a £(3)-dependence that cancels in the sum. 
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of diagram (/) in the last section. Including the respective symmetry factors ^, ^, | we 

get «-(i) 5 M^ in3 (a-I'<g^K9-3' (70) 

For future use we calculate the improper four loop <^ 4 -diagram (t) where the dotted line 
means a '(— l)-fold' propagator (x — y) +2 . The result is (with a symmetry factor of ^) 

«-(f) 4 p(^(S)-i"(S) + i ta (S)-S^«)- < 73 > 



Integration by parts. We determine the following symmetry factors: (m): ~, (n): ~, 

(o): |, (p): |, (?): §, (r): ±. 

The idea is to use Eq. fl65|) to relate the above graphs among each other. Sometimes 
it will be necessary to multiply Eq. (|65|) by (x\ — x 2 ) 2 , (x\ — x 3 ) 2 , or (x 2 — x 3 ) 2 . Since 
these factors are independent of x they do not affect partial integration with respect to x. 
However, if these factors do not combine with propagators (x\ — x 2 )~ 2 , etc., one obtains 
improper y9 4 -graphs like diagram (£). In most cases it is possible to eliminate those graphs 
by a second application of Eq. (|65|). In the following table we denote first the graph we 
start from, then the variables which correspond to (x, X\, x 2 , x 3 ) in Eq. (|65D (according 
to Fig. 1), occasionally the variables of a second application of Eq. (|65|), and finally the 
resulting equation including the symmetry factors. 

graph (x, xi, x 2 , x 3 ) (x,xi,x 2 ,x 3 ) equation 
~~ (p) (2,1,4,3) = G p + G = 2G m - ±gG g 

(r) (3,4,1,2) \G r = G q -\gG h 

(r) (3,4,2,1) (2,4,1,5) \G r + 2G n + 2G m = 4G dA - \gG h - gG d , 3 

(m) (2,3,5,1) (4,3,1,5) 2G q + G Q + 2G m = 4G dA - \gG h - gG<u, 

~~ 5) WAZ) - Gh = -G g + 2G d;3 - \gG d ,2 (Eq. ©) 

G?) (2,1,3,4) G g = 4G t - \gG d , 2 

gG t = 12G k + ±G d , 2 (Eq. ©, (0), ©) 

(74) 

The last equation can explicitly be checked by looking at the amplitudes. We recognize 
that there are only four equations to evaluate six five loop diagrams. However summing 
up the first four equations gives 

2G m + 2G n + 2G + G p + G q + G r = 8G dA - (g/4) ■ (2G g + 5G h + 8G d>3 ) . (75) 

With the last three equations in the table we can express the left hand side completely in 
terms of Fourier amplitudes of the y2 4 -theory 

2G m + 2G n + 2G + G p + G q + G r = 8G dA + 3QG k - 18 (g/4) G d , 3 + y (g/4) 2 G d>2 
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(g\ 5 l(l 27 3 /p 2 \ 299 2 /p 2 \ 809 ( p 2 \ 1853 \ 

(76) 

The graphs (m), (n), (o) are not symmetric under interchange of the external legs. There- 
fore we have to count them twice in the two-point function and the left hand side becomes 
exactly the combination we want to calculate. 

The bipyramide graph. The bipyramide graph (s) is the first non-planar two-point 
graph and commonly regarded as the most complicated five-loop diagram. It was first 
calculated in 1981 within dimensional regularization by K.G. Chetyrkin and F.V. Tkachov 
Recently it was analyzed within differential renormalization by V.A. Smirnov [IT]. 



So it is a good candidate to test the power of our calculation scheme. We work in 
coordinate space. It is convenient to introduce a quaternionic notation. The inversion of 
a quaternion x is given by x \— > - which can be understood in the four vector language as 
inversion of the length of x (\x\ i— > |x| _1 ) and a reflection at the z-axis (the direction of 
the unit quaternion 1). The square (x — y) 2 becomes the square of the absolute \x — y\ 2 , 
however we stick to the brackets in the following calculation to keep the notation more 
transparent. 

The variables (x,a,b,c,y) correspond to (1,2,3,4,5) in Fig. 1. We have to calculate 
the following integral 

dci/Cl/bcl/C 2 2 2 2 2 2 2 2 2' 

(x — a) (x — b) [x — c) (a — b) (a — c) (6 — c) (a — y) {b — y) (c — y) 

(77) 

The external legs are amputated, they can easily be added in the end. 

The integral is convergent at infinity (it is logarithmically divergent at a = b = c = y 
and a = b = c = z) and therefore the integration variables a, b, c can be shifted by y. 
With z = x — y we have 

/■ll 1 

I (z) = / dadbdc—T—— ^ ^ t ^ ^ s- . (78) 

v ' J aW& (a - b) z (a — c) (b — c) (z — a) (z — b) (z — c) 

With the inversions a' = b' = |, c' = ~ one obtains (d 4 a = a'~ 8 d 4 a') 
da'db'dc' 



a'%' 6 c m 



( Z a') \ Z If) \ Z c') (a' b') (a' c') (v c> ) 



^ da'db'dc' 



2 • 



(a' - ±y (b' - I)' (c' - I) Z " a'Y (C - a') z (C - 6') 
A shift a" = a' - |, 6" = &'-£, c" = c' - ± yields 

4 / da"db"dc" J- 5 5 r . (79) 

z eJ a " 2 b" 2 c" 2 (a" -b"f(a" -c"f(b" -c"f V ' 
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It seems that we have lost the z-dependence in the integral. However, since the integral 
is still divergent, this is not the case as we will see soon. 
We finally use the rescaling b" = a"u, c" = a"v to obtain 



da" 
a" 1 



du I dv 



v) 2 v 2 



(1 — u) v? (u 

The u- and f-integral is finite and gives a positive number 



(80) 



Gegenbauer polynomial techniques (e.g. ||12 |, 



It can be evaluated using 
with the result! |C(3). The angular 
integral in / da" /a" 4 gives 2n 2 . Including the normalization (Eq. (|35|)) one is left with 
the radial integral \ d\a"\/\a"\ = \ 1h(|o // |)|°//^q°. If we would have started with an 
integral in the a"-variable this integral would give zero. However as discussed at the end 
of Sec. |2.2j it is now essential to reintroduce the original variable a before approaching the 
limits. Since a" — - — 



Jd\a"\ 


-H 


i 


1 


) 


J 2\a"\ 




a 


z 




o 1 1 









a=0 




z — a\ 
lazl 



a=0 




Collecting all pieces we have finally found 

I(z)= 3 




(81) 



(82) 



The transformation to momentum space is given by Eq. (0). Including the external legs 
and the symmetry factor | we obtain 



G r . 



P 



P 
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7 3UC(3)ln a "SCO) In + (3) 



P 



A 2 



A 2 



4 



33) 



Comparing with dimensional regularization |TIJ gives the minimum coincidence that both 
results are proportional to C(3)- It is not possible to be more precise since in 0] only the 
singular part was calculated. Note that the techniques we used can not be generalized to 
dimensions different from four. 

It is also hard to compare our result with the one gained by differential renormalization 
in since the author restricted himself to regularize the amplitude and did not evaluate 



the rather complicated integrals over the internal variables. 
3.4.5 The two-point Green's function 

Collecting all results from the last sections we obtain for the full propagator of the y? 4 - 
theory 



G(P) = £ (l + (f 




HIM 



6 Most efficiently one uses the identity \uv\/(u — v) 2 = — dP (\u\ / \v\) lP Yl^Li nC n -i(u ■ v) ■ {n 2 + 
P 2 )^ 1 and the orthogonality relation — J dxnC n -i(y • x)mC m -\(i ■ z) = S n m nC n -i(y ■ z) to obtain 
iroc d pv°° 1 n 2 (n 2 +F 2 )- 3 . 
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9V ( 3 ,..3 [P 2 \ 59 2/V\ , 1535, /V\ 121 \ 



A) \T6 ln {Ti)-36 ln ^J+W^^J-lsJ 




4/ \32 \A 2 / 288 \ A2 
Comparison with differential renormalization || 

^>)^( 1 + (f)^ ln (x) 2+ (f) 3 (i b2 (i) 2 -i ln (f)> -) w 

shows that only the leading logarithms coincide. Notice that one never gets ln-independent 
terms in 0. 

If A is rescaled according to In A i— > In A + 1 in Eq. (|85|) the logarithmic terms coin- 
cide with that of Eq. (153). The ln-independent terms can be adjusted via a momentum 



independent rescaling by 1 — j^(g/A) 2 + j^(g/A) 3 . However, at this point it is not clear 
whether the differences disappear after appropriate redefinitions also at higher orders. 

3.4.6 The renormalization group 

It is possible to extract the /3-function and the anomalous dimension 7 from the two-point 
function alone if one assumes that j3 and 7 are independent of A. Moreover the existence 
of a renormalization group equation is a non-trivial test for the renormalization scheme. 
Comparing the coefficients in 

' 9 +?f!]i + ^!))cf?,A, P ] = » (86) 



yields 



<91nA r \AJ d(g/A) \AJ I \A 



"(f) - 3 (f) 2 -?(!) 3+ (? +12 « 3 »)(^ +0 (^)' < 87 > 



A) 12 \AJ 96 \AJ 192 V4 

The first two terms of f3 and the first term of 7 are standard. The coefficient in front of the 
C(3)-term also coincides with other schemes [0, [ 13| . However e.g. the vanishing third or- 



der and the £(3)-independent fifth order term of 7 is specific to our scheme. In differential 
renormalization § one obtains f3(g/A) = 3(g/A) 2 - f {g/Af + (31 + 12C(3))(#/4) 4 + . . ., 
7(5/4) = ^/4) 2 -§(s/4) 3 + .... 



4 Results and outlook 



A new renormalization scheme was proposed. It provides all amplitudes fully renormal- 
ized, it has no explicit cutoff or counterterms and allows to keep the spacetime dimension 
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fixed. The scheme defines all integrals in an unambiguous way, it thus corresponds to a 
definite choice of a subtraction prescription. 

The renormalization scheme emerges from differential renormalization by an a priori 
fixing of all integration constants at their mathematically most natural values. It is closely 
related to the theory of generalized functions. 

We demonstrated how to use this renormalization scheme if applied to the toy prob- 
lem of a two-point (mass) interaction in coordinate space. Although this theory is non- 
renormalizable by power-counting it was possible to recover the correct result within our 
scheme. A theorem was presented that allowed us in a more general framework to recon- 
struct the full result from such a singular expansion. With this theorem it was possible 
(but more complicated) to regain the true result even for the dimensionally regularized 
toy model which failed to give the correct perturbation series. 

The main application of our scheme was the y9 4 -theory. Equations that are very special to 
four dimensions and to our renormalization prescription enabled us to calculate the two- 
point Green's function up to five loops (Eq. (|34"D). Most remarkable was the observation 
that at (four) five loops the diagrams are organized in such a way that a (one-) two- 
fold underdetermined system of linear equations could be solved for the sum over certain 
diagrams. This made the evaluation of many single graphs needless. 

We were left with the nonplanar five-loop graph which could as well be calculated 
analytically in our renormalization scheme. It is obvious that the dimension of spacetime 
plays a crucial role in the calculation of the bipyramide graph (as it does for the matching 
of diagrams via integration by parts). Only in four dimensions the coupling becomes 
dimensionless. The resulting conformal symmetry was used via the inversion a h- >• a -1 as 
the most essential step in the evaluation of the integrals. 

The two-point function was compatible with the renormalization group and it was 
possible to extract the /3-function up to fourth and the anomalous dimension 7 up to fifth 
order in the coupling. 

For future work the idea of grouping certain classes of diagrams and calculating their 
sum without referring to single graphs appears especially promising to us. We expect 
that the matching of amplitudes persists to some extent at higher orders. In this way 
perturbation theory could be simplified and even analytical results beyond the fifth order 
may be possible. (Recent calculations confirm this for the sixth order of y9 4 -theory.) Most 
desirable would be to find the general structure that organizes the amplitudes to groups 
that can be evaluated via integration by parts. General questions of renormalizability and 
the problems related to the multiplication of generalized functions have to be investigated 
more carefully. 

A goal of obvious importance is the application of natural renormalization to gauge 
theories. In general one has to avoid conflicts between Ward identities (reflecting gauge 
symmetry) and the renormalization scheme. This problem is already present in two 
dimensions and can be solved by using the transverse (Landau) gauge. The Schwinger 
model can be solved within this framework by summing up the whole perturbation series 
(e.g. the fermion correlation function) ||. The key tool is, similar to the </2 4 -theory, to 
calculate whole classes of Feynman diagrams without evaluating single amplitudes. In 



REFERENCES 



23 



four dimensions first results are promising, however for QED we have not yet found how 
to group diagrams to simplify the calculations. 
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